Draft version February 5, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. ^jllj^A, 



QUASARS AT Z=6: THE SURVIVAL OF THE FITTEST 
Marta Volonteri^ & Martin J. Rees^ 



o 
o 

(N 



in 



> 
m 
a^ 
o 
l> 
o 
\o 
o 



Oh' 

I 

o 



X 



(Received 



Accepted 



Draft version Pebruary 5, 2008 

ABSTRACT 

The Sloan Digital Sky survey detected luminous quasars at very high redshift, z > 6. Follow-up 
observations indicated that at least some of these quasars are powered by supermassive black holes 
(SMBHs) with masses in excess of 10^ M0. SMBHs, therefore, seem to have already existed when 
the Universe was less than 1 Gyr old, and the bulk of galaxy formation still has to take place. We 
investigate in this paper to which extent accretion and dynamical processes influence the early growth 
of SMBHs. We assess the impact of (i) black hole mergers, (ii) the influence of the merging efficiency 
and (iii) the negative contribution due to dynamical effects which can kick black holes out of their 
host halos (gravitational recoil). We find that if accretion is always limited by the Eddington rate 
via a thin disc, the maximum radiative efficiency (spin) allowed to reproduce the LF at z = 6 is 
e = 0.12 (a = 0.8), when the adverse effect of the gravitational recoil are taken into consideration. 
Dynamical effect unquestionably cannot be neglected in studies of high-redshift SMBHs. If black 
holes can accrete at super-critical rate during an early phase, reproducing the observed SMBH mass 
values is not an issue, even in the case that the recoil velocity is in the upper limits range, as the 
mass ratios of merging binaries are skewed towards low values, where the gravitational recoil effect is 
very mild. We propose that SMBH growth at early times is very selective, and efficient only for black 
holes hosted in high density peak halos. 

Subject headings: cosmology: theory - black holes - galaxies: evolution - quasars: general 



1. INTRODUCTION 

It seems a challenge for theoretical models to explain 
the luminosity function of luminous quasars at z » 6 i n 
the Sloan Digital Sky Survey CSDSS iFa.n eraTI 12001 Ii . 
The luminosities of these quasars, well in excess of 10 
erg/s, imply SMBHs with masses IQ^Mq already in 
place when the Universe is only 1 Gyr old. The accretion 
of mass at the Eddington rate causes the black hole (BH) 
mass to increase in time as 



M{t) = M(0) exp 



1 - e t 



(1) 



where isdd = 0.45 Gyr and e is the radiative efficiency. 
Among the seed BHs proposed, the more commonly in- 
voked (e.g. PopIII star remnants, gravitationally col- 
lapsed star clusters) have masses in the range M(0) = 
10^ - lO^M©, forming at z = 30 or less. Given M(0), 
the higher the efficiency, the longer it takes for the BH to 
grow m mass by (say) 20 e-foldings l)ShaDird !2005'). For 
a Schwarzschild black hole, the standard thin disc radia- 
tive efficiency is e w 0.06, and there is plenty of time for 
the BH seed to become supermassive. The timescale to 
grow from M(0) = 10^ - lO^M© to ~ lO^Mg is less than 
0.5 Gyr. 

If accretion is via a geometrically thin disc, though, the 
alignment of a SMBH with the angular momentum of the 
accretion disc tends to efficiently spin holes up (Volon- 
teri et al 2005), and radiative efficiencies can therefore 
approach 30-40%. With such a high efficiency, e = 0.3, 
it takes longer than 2 Gyr for the seed to grow up to a 
billion solar masses. 

It seems therefore difficult to reproduce the observa- 
tional constraints without invoking exotic processes. 
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We take here a conservative approach and try to 
determine the parameter space which allows a SMBH 
growth compatible with observational constraints at z ~ 
6, as evinced from the luminosity function of quasars 
l|Fan et al]l2004l) and observations of galaxies at z ~ 6. 
In this paper we critically assess models for the early evo- 
lution of SMBHS, exploring the parameter space of the 
involved processes: accretion rate, radiative efficiency, 
dynamical processes and the initial density of massive 
black hole (MBH) seeds. 

At z < 5 MBH mergers do not play a fundamental role 
in building up the mass of SMBHs (Yu & Tremaine 2003), 
but they can be possibly important at z > 5, where we do 
not have constraints from a Soltan-type argument, which 
compares the local MBH mass density with the mass 
density inferred from luminous quasars, as the luminosity 
function of quasars is not constrained at z > 6. Mergers 
can possibly contribute positively to the build-up of the 
high redshift SMBH population (Yoo & Miralda-Escude 
2004), as they contribute to make a big black hole from 
many small seeds. 

On the other hand, dynamical processes can disturb 
the growth of BHs, especially at high redshift (Haiman 
2004, Yoo & Miralda-Escude 2004), and give a negative 
contribution to SMBH growth. In the shallow potential 
wells of mini-halos, the growth of MBHs can be halted 
by the 'gravitational rocket' effect, the recoil due to the 
non-zero net linear momentum carried away by gravita- 
tional waves in the coalescence of two black holes. Also, 
if MBHs are widespread, and binary black holes coales- 
cence timescales are long enough for a third MBH to fall 
in and interact with the central binary, Newtonian three- 
body interactions can lead to the expulsion, or recoil, of 
the binary. The accretion history must then be studied 
jointly with the dynamics involving MBHs. 
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Yoo & Miralda-Escude (2004) explored the minimum 
conditions that would allow the growth of seed MBHs 
up to the limits imposed by the highest redshift quasar 
currently known: SDSS 1148+3251. This quasar, at 
z = 6.4, has estimates of the SMBH mass in th e rang e 
(2 - 6) X IO^Mq lIBaTth et a,1. 2003; Willott et a,1Jl200l . 
Yoo & Miralda-Escude (2004) showed that the mass 
of this SMBH can be explained assuming (i) continued 
Eddington-limited accretion onto MBHs forming in halos 
with Tvir > 2000K at z < 40. Their model assumes, also, 
(ii) a high influence of BH mergers in increasing the MBH 
mass: a contribution by itself of order 10^ M0. Their 
investigation takes into account the negative feedback 
due to the 'gravitational rocket' effect (see also Haiman 
2004). Recent estimates suggest a more modest recoil 
velocity, compared to the typical values adopted in the 
past. We test here how much difference (if any) these 
new estimates imply for the growth of black holes in pre- 
galactic halos. 

It is also important to understand where z = 6 quasars 
are hosted. By matching the number density in haloes 
more massive than Mh. to the space density of quasars at 
z = 6 l|Fan et al.ll2fffll . Mh = IO^^Mq. This assump- 
tion corresponds to requiring that the duty cycle of high 
redshift quasars is unity, i.e. all the BHs inhabiting halos 
with mass larger than 10^'^ Mq are a ctive. The only avail- 
able observation of a quasar host ijWalter et al. 2004) , 
shows interesting features. The kinematics of the ob- 
served molecular gas implies the lack of a massive bulge 
around the SMBH, but suggests a DM halo with a mass 
similar to that of the predicted bulge, of order 10^^ M0. 
This is much less than = 10^'^ Mq and allows for a 
much smaller duty cycle. The uncertainties on the dy- 
namical configuration of the gas (e.g. inclination) are 
still large, however. 

In the next section we review the model for assembly 
of MBHs in cold dark matter (CDM) cosmogonies. We 
then review models for MBH growth by accretion in high 
redshift halos (§3), and discuss the dynamical evolution 
of MBH binaries in pre-galactic systems (§4). The results 
from the interplay of accretion and dynamical processes 
are presented in §5. Finally, we discuss the implications 
on the global evolution of the SMBH population (§6). All 
results shown below refer to a ACDM world model with 
= 0.3, f^A = 0.7, h = 0.7, f7b = 0.045, erg = 0.93, 
and n = 1. 

2. ASSEMBLY OF PREGALACTIC MBHs 

The main features of a plausible scenario for the hier- 
archical assembly, growth, and dynamics of MBH s in a 
ACDM cosmo logy have b een discussed in Volontcri et al.l 
l|2003l f200,^ : iMadau et all l|2004D . Dark matter halos 
and their associated galaxies undergo many mergers as 
mass is assembled from high redshift to the present. The 
halo merger history is tracked backwards in time with 
a Monte Carlo algorithm based on the extended Press- 
Schechter formalism. "Seed" holes are assumed to form 
with intermediate masses in the rare high v — a peaks 
collapsing at z = 20 - 25 ()Madau & R,ee^l200H) as end- 
product of the very first generation of stars. 

As our fiducial model we take // = 4 at z = 
24, corresponding to a mass density parameter in 
MBHs of order p, ~ lO^MoMpc^^, about three or- 
ders of magnitude smaller than the local SMBH mass 



densitv JXller fc Richstol^ 120021: iMerloni et all l2004t 
iMarconi et all 120041 lYu fc Tremaindl2002D . The mass 
density in seed MBHs cannot be larger than the mass 
density of SMBHs at z = 3, pbh{z = 3) ~ 4.5 x 
lO'^MoMpc-^ (Merloni 2004). 

This choice of se ed MBHs occupation is similar to that 
of iVolonteri et aP |2003) (seed holes in 3.5cr peaks at 
z = 20) , ensuring that galaxies hosted in halos with mass 
larger than 10^^ M© are seeded with a MBH. The as- 
sumed threshold allows efficient formation of SMBHs in 
the range of halo masses effectively probed by dynamical 
studies of SMBH hosts in the local universe. We will also 
consider a case in which seed holes are more numerous 
(e.g. 25 times more than the fiducial case). 

We are here interested in the evolution of the uttermost 
massive halos pres ent at z = 6, i.e. very highly biased 
structures (Dicman d et al.ll20 05). As a consequence, the 
density of seed holes within the volume occupied by the 
progenitors of the halo is large enough that the merging 
of two minihalos both hosting a BH is not a rare event. 
This is in contrast with the average density of seed holes. 
In a cosmic volume, the seed h oles typically evolve in 
isolation fcfr. IMadau et al.ll200l) . 

We generate Monte Carlo realizations (based on the 
extended Press-Schechter formalism) of the merger hi- 
erarchy of a M;i = 10" Mq halo at z = 6. The halo 
mass is chosen by requiring that the number density in 
haloes more massive than Mh match es the space density 
of quasars at z = 6 ijFan et al .1120041) . 

3. ACCRETION AND RADIATIVE EFFICIENCY 

We explore here two regimes: in one case the ac- 
cretion rate is limited to the Eddington rate, in the 
second case MBHs are allowed to accrete at super- 
critical rates (cfr. IVolonteri fc R,ee3 l|2005l) for a de- 
tailed description of the model). In all cases we assume 
that accretion episo des ar e triggered by major rnergers 
iMihos & Hernauist'lggl, TflOffl iDi Matteo et al.ll2005t 
^azantzidis ct al. 2005), which we define as mergers be- 
tween halos with a mass ratio 1:10 or higher. Mergers 
with smaller mass rati os are pro bably unable to trigger 
substantial gas inflow l)Coxll200^ . 

In the first case, the hole accretes at the Eddington 
rate a gas rest mass Amg. This leads to a change in the 
total mass-energy of the hole given by 

Am = 7.7 X 10^ Mq V^foo^ (2) 

where Vc,ioo is the circular velocity of the merged sys- 
tem in units of 100 kms~^. Adopting Equation El im- 
plies assuming that the correlations between black hole 
mass, velocity dispersion and circular velocity arc main- 
tained th roughout cosmic time ( Tre maine et al. 2002|: 
lFerrares'3 ^002') . The quantities Am and Atoq are re- 
lated by Am — (1 — e) Amp, where e is the mass-to- 
energy conversion efficiency, equal for thin-disk accre- 
tion to the binding energy per unit mass of a particle 
in the last stable circular orbit^. The MBHs spin, S, 
is modified dur i ng th e accretion phase as described in 
IVolonteri etldl We are interested here in deter- 

The simple relation would be modified when the thickness of 
the disc is of order its radius and can also changed by magnetic 
effects whi ch allow energy rele ase from within the innermost stable 
orbit fe.g.. IKrolik et alil2005^ 
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mining the maximum radiative efficiency (or spin, us- 
ing tlie standard conversion) which allows the growth of 
MBHs into SMBHs with masses larger than 10^ Mq be- 
fore z = 6. As a reference, we have considered an upper 
limit to the radiative efficiency of e = 0.16. This corre- 
sponds, adopting the standard conversion for accretion 
from a thin disc, to a maximum spin parameter of the 
BH g = 0.9. This v alue was chosen in agreement with 
IGammie et al.l l^2004^ simulations, which suggest that the 
maximum spin MBHs can achieve by coupling with discs 
in magncto-hydrodynamical simulations is a ~ 0.9. 

In the second case, we focus only on a subset of halos, 
those with effective atomic cooling. The cooling curve of 
metal-free gas has a sharp break at T < IC^K, so that 
for halos with Tvir < lO^K the only available coolant 
is molecular hydroge n. is nevertheless a rather in- 
efficient coolant iMadau et alJl2001l) . Metal-free halos 
with virial temperatures Tvir > 10 K can instead cool 
even in the absence of H2 vi a neutral hydrogen at omic 
Hues to - 8000 K. Following lOh fc HaimaiJ l|2002t) . we 
assume that a fraction of the gas settles into an 
isothermal, exponent ial disc, embed d ed in a dark matter 
halo described by a iNavarro et al.l l)1997|) density pro- 
file. The mass of the disc can therefore be expressed as 
Mdisc — fd{^b/^M)Mh- The disc is geometrically thick 
and i t has a very high central de nsity (c fr Bro mna fc Lc 
120031) . We refer the reader to^^^teri fc ReST pC 
and i Lodato fc NataraianI l)2006j) for a detailed descrip- 
tion of the model. 

We estimate the mass accreted by the MBH from 
the surroun ding disc within the Bondi-Hoyle formalism 
lIBondi fc Hovlelll94l . 



Bondi " 



a 47r ^miin 



= 4 X 10" 



M, 



BH 



10^ Mr. 



V104cm-3/ V 



8000K 



-1.5 



M, 



(3) 



-3 < 



where uq is the density of the gas, of order 10"^ cm 
no < lO^cni^'^ at the center of the gas discs in high 
redshift halos. 

The collapsing gas disc likely rotates as a rigid body, 
rotation is therefore small during the initial collapse 
phase, and the infall of gas is quasi radial. The size of 
the accretion disc, rjn, is of order of the trapping radius: 



M 



rtr 



rs 



M\ 



BH' 



(4) 



Edd 



i.e. the radius at which radiation is trapped as the in- 
fall speed of t he gas is lar ger than th e diffusion speed of 
the ra diation. iBegelmanI l|1979j) and iBegelman fc Meieil 
1)1982(1 studied super-critical accretion onto a BH in 
spherical geometry and quiescent thick discs respectively. 
Despite the uncertainties, it still seems possible that 
when the inflow rate is super-critical, the radiative ef- 
ficiency drops so that the hole can accept the mate- 
rial without greatly exceeding the Eddington luminosity. 
The efficiency could be low either because most radia- 
tion is trapped and advected inward, or because the flow 
adjusts so that the materi al can plunge in from an orbi t 
with small binding energy l|Abramowicz fc Lasotall980|) . 

The accretion rate is initially super-critical by a factor 
of 10 and grows up to a factor of about 10* (Volonteri fc 



Rees 2005), thus making the flow more and more spheri- 
cal (see Equation^ . On the other hand, the radius of the 
accretion disc increases steeply with the hole mass, until 
the whole 'plump' accretion disc grows in size enough to 
cross the trapping surface. We assume here that, when 
the radius of the accretion disc becomes a factor of 5 
larger (although the choice of the exact value is somehow 
arbitrary) an outflow develops, blowing away the disc. 
Typically, assuming this threshold, a MBH accretes only 
a fraction (1-10%) of the mass in the disc. Also in this 
ca se we follow the evolution of MBH spins as described 
in lVolonteri~ al. (2005). 

An alternative model for super-critica l accre tion on 
seed MBHs can be found in iBegclma n et alJ (|2006) , 
where unstable discs are considered, instead of stable sys- 
tems. Both models for supercritical accretion apply only 
to metal-free halos with virial temperature Tvir > 10'*K, 
that is, to rare massive halos, at the cosmic time consid- 
ered. Rapid early growth, therefore, can happen only 
for a tiny fraction of MBH seeds, in a selective and 
biased way. For illustrative purposes. Figure ^ shows 
how different accretion rates, or the same accretion rate 
with different radiative efflciencies can dramatically mod- 
ify the mass of a MBH starting accreting a,t z — 24. 




Fig. 1. — Growth of a MBH mass under different assumption 
for the accretion rate and efficiency. Eddington limited accretion: 
e = 0.1 {solid line), e = 0.2 (short dashed line), e = 0.4 (dot-dashed 
line). Super-critical accretion, as in Volonteri & Rees 2005 (long 
dashed line). 



4. DYNAMICAL EVOLUTION OF MBHs 

To include in this study the dynamics of MBHs, we 
will follow the evolution of MBHs and t heir hosts in full 
detail with a semi-analytical technique l|Volonteri et alJ 
I2OO3. 20051. 

4.1. MBH binaries merging efficiency 

The merging - driven by dynamical friction against 
the dark matter - of two comparable-mass halo+MBH 
systems ( "major mergers" ) drags in the satellite hole to- 
wards the center of the more massive progenitor, lead- 
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ing to the formation of a bound MBH binary with sep- 
aration of ~ pc. In massive galaxies at low redshift, 
the subsequent evolution of the binary may be largely 
deter mined by three-body in teractions with background 
stars ijBegelman et al.l ll98rf). Dark matter particles will 
be ejected by decaying binaries in the same way as the 
stars. Another possibility is that gas processes, rather 
than three-body interactions with stars or DM, may in- 
duce MBH b i naries t o shrink rap idl y an d coalesce (e.g. 
iMaver ct al. (2006^; 'Dotti ct al. (2006); Escala et alj 
12004): Li u ct al. (2003); Armitagc & Nataraian (200% 
iGould fc R ix ( 2000)). If stellar dynamical or gaseous 
processes drive the binary sufficiently close (< 0.01 pc), 
gravitational radiation will eventually dominate angular 
momentum and energy losses and cause the two MBHs 
to coalesce. 

In gas rich high redshift halos, the orbital evolution of 
the central SMBH is likely dominated by dynamical fric- 
tion against the su rrounding gaseous medium. The avail- 
able si mulations f Maver et all I2006t IDotti et all I2006t 
lEscala ct al. 2004) show that the binary can shrink to 
about parsec or slightly subparsec scale by dynamical 
friction against the gas, depending on the gas thermo- 
dynamics. These binary separations are still too large 
for the binary to coalesce within the Hubble time owing 
to emission of gravitational waves. On the other hand, 
the interaction between a binary and an accretion disc 
can lead to a very efficient transport of angular momen- 
tum, and the secondary MBH can reach the very sub- 
parsec separations at which emission of gravitational ra- 
diation dominates on s hort timescales Liu et al. ( 2003) 
lArmitaee &: NataraianI 



ifMol : ICxOuld fc Rixl ll200r 

The viscous timescale depends on the properties of the 
accretion disc and of the binary: 



ivis =0.1 Gyra^^^ 



-1 / 



mi 



V lO^Mc 



-1/2 



(5) 



where ai is the initial separation of the binary when the 
secondary MBH starts interacting with the accretion disc 
in units of parsec, (h/r) is the aspect ratio of the accre- 
tion disc, h/r = 0.1 above, a is the Shakura & Sunyaev 
viscosity parameter, a = 0.1 above, and mi is the mass 
of the primary MBH, in solar masses. The emission of 
gravitational waves takes over the viscous timescales at 
a separation ijArmitage fc Nataraianll2?)(i5l) : 



aaw = 10 f- 



-16/5 



0.1 



-8/5 3/5 / 

"o.i qo.i 



nil 
^lO^M, 



-0 . 



(6) 



where q = 1712/ mi < 1, is the binary mass-ratio. The 
timescale for coalescence by emission of gravitational 
waves from acw is much shorter than the Hubble time: 



bc>a\t) 



256G^mim2{mi + 7712) 



(7) 



We have assumed here that, if an accretion disc is 
surrounding a hard MBH binary, it coalesces instanta- 
neously owing to interaction with the gas disc. If instead 
there is no gas readily available, the binary will be losing 
orbital energy to the stars, if the initial mass function 
of PopIII stars is bimodal a nd a large pop ulation of low 
mass stars is formed (,,Nakamura fc Umemura 2001) or 
to dark matter background otherwise. The coalescence 



timescale typically becomes m uch longer than the Hubble 
time (see IMadau et al.l (j2004j) for a thorough discussion 
for this case). 

To test the influence of the merging efficiency, we have 
compared the above merging efficient model to a conser- 
vative merging inefficient case in which the interaction 
with gas is neglected, and the binaries shrink only via 
DM-|-stars scattering, under the assumption that the loss 
cone stays full. 

4.2. Gravitational recoil 

At high redshift, the recoil as a result of the non- 
zero net linear momentum carried away by gravitational 
waves, may cause the ejection of M BHs from the shal- 
low potential wells of their hosts (e.g. .IMadau e t al. 200^ 
IMadau fc Quataerti2004l:lMerritt et alJ2004j) . The recoil 
velocity has still large u ncertainties. E arly calculations 
in the Newtonian regime ljFitchettlll983|) . predict at most 
a recoil velocity of w lOOkms"^. The Newtonian calcu- 
lations, in the circular case, predict the center of mass 
to move on circular orbits, spiraling outward while the 
binary orbit spirals inward, with a velocity 



VCM = 1480 km s 



fil) 



fn 



Rl 



(8) 



where Rs = 2G{mi+m2)/c^ is the Schwarzschild radius 
of the system, and the scaling function reads f{q)/fma^, 
with fiq) = g2(l _ g)/(l + and /„,ax = 0.01789. 

Here Rl represents the radius of the last stable circular 
orbit, which is Rl = QGmi/c^ — ?>Rsmi/ (mi -I- 7712) in 
the Schwarzschild metric. 

The recoil velocity during the plunge phase probably 
has the largest contribute. Calculations of the gravita- 
tional recoil inside the innermost stable orbit (ISCO) 
naturally have large uncertaintes, but the recoil veloc- 
ity should be constrained between the upper- and lower- 
limits (Vupper and Viower respectively) suggested by Fa- 
vata et al. (2004), which span a range between a few 
kms"""^ and several hundreds kms~^ for a binary with a 
mass ratio q cii 0. 1. 

iBlanchet et al.l 1)2005') calculated the gravitational re- 
coil at the second post-Newtonian order for non-spinning 
holes. Their calculations, available at the moment 
for Schwarzschild holes only, narrow the uncertain- 
ties on the final veloc ity from more than an order of 
magnitude llFavata et al . 2004), down to about 50%. 
iDamour fc Gopakumail ((2006.) . using the effective one- 
body approach, predict velocities about a factor of 3 
less than Blanchet et al. The uncertainties quoted 
in Blanchet ct al. (200^ and fDamour fc GoDakuma'3 
(j2006) make the results of the two calculations incom- 
patible. 

The latest estimate on the rec oil comes from fully rel- 
ativistic numerical simulations ({Baker et al.l 1200(11) fol- 
lowing the dynamical evolution of a black hole bi- 
nary within the ISCO. These simulations, carried for 
a m ass ratio q = . 2, pre dict a recoil midway through 
the IBlanchet et aP l)2005D and IDamour fc Gooakumai 
f2006') predictions. The recoil predicted bv lBaker et al. 
(2006) is still large enough to eject the merging bi- 
nary from small pre-galactic structures. As shown by 
Yoo fc Miralda-Escude 2004 , if th e recoil effect is mild 
(e.g. IDamour fc Gopakumail l)2006j) : lFavata et all l|2004fl 
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simistic (although reahstic) conditions under which it is 
stiU possible to develop the population of z = 6 quasars, 
we adopt as a zeropoint the Ba ker et alJ l)200 6l results. 
We rescale the recoil to different mass ratios adopt- 
ing Fitchett's scaling function. Although Fitchett de- 
rived the scaling in the perturbative regime, outsid e the 
strong gravity region, both Blanch et et alJ lj2005r and 
iDamour fc Gopakumail l|2006|l find a dependence on the 
mass ratio in very good agreement with Fitchett's scal- 
ing. 



!,1000 F- 



lower limits), it is easy to fulfill the constraints on quasars altogether. If instead dcm < I'osc, the binary re- 
at z = 6'^. As we are considering here the most pes- coils within the galaxy and it s orbit will (slowly) decay 

owing to dynamical fric tion ijMadau fc Quataerj l2004t 
IVolonteri fc Perna j |200,^ . The escape velocity calcula- 
tion assumes an NFW halo for the DM component. The 
baryonic component is modeled as an isothermal sphere 
truncated at the radius of the MBH sphere of influence, 
except in the cases where gas can condense in a disc (see 
section 3). In this case the escape velocity is computed 
assuming an exponential profile for the gaseous disc. The 
escape velocity is between 2.5 — 5 times the virial velocity 
of the host halo. Note, however, that if the recoil veloc- 
ity is larger than about twice the velocity dispersion of 
the host halo, the dynamical friction timescale for the 
ejected MBH to return to the center of the halo is larger 
than the Hubble time (Madau fc Quataert 2004). 

5. RESULTS 

The combination of a halo merger tree and our semi- 
analytical scheme to trt the growth of MBHs and their 
dynamics is a powerful tool for tracking the evolution of 
MBHs at early times. Our semi-analytical approach is 
highly idealized, but allows us to explore a large range 
of different scenarios and their consequences on the evo- 
lution of the MBH population at high redshift. 

5.1. MBH dynamics 

The dynamical and accretion evolution of high-redshift 
MBHs is so strongly intermingled that it is very difficult 
to clearly separate the relative effects. The accretion 
history determines the MBH binaries mass-ratios, and 
therefore the merging timescale, and the recoil velocity 
(cfr. eqEl). 
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Fig. 2. — Merging BHs recoil velocity as a function of bi- 
nary mass-ratio. Solid line: spin parameter a = 1. Long-dashed 
line: spin parameter a = . D ot-dashed line: spin parameter 
a = —1. Bottom-left: IFavata e t al. 12004) upper limits. Top- 
left: IBalc er et al.l 120061) . The mass ratio dependence is estimated 
by the Fitchett's formula, and the spin-orbit i nterac tion as in 
ll'avata et al. 1.2004f). Top-right: IBlanchet eTaD 120051) . Bottom- 
riaht: IPamourT^jOpakumarl |200^ 



One remaining issue is the effect of black hole spins, 
and the spin-orbit coupling^ which however, seems to 
be mild. IFavata et alJ lj2004f) suggest a modification of 
Fitchett's scaling function to keep into account spin- 
orbit coupling: f{q) = f{q)\l + (7/29)ai/(l - q)\/\l + 
(7/29)di/(l — 0.127)1. Here cii is the physical spin pa- 
rameter of the large hole (5*1 — SiGMbh/c, < oi < 1), 
and a is the "effective" spin parameter of the binary sys- 
tem. Damour suggests a = ai(l + 3g/4)/(l -I- g)^ in the 
post- Newtonian limit. Figure |2l compares different theo- 
retical estimates for the recoil velocity. 

We therefore adopt the following expression to param- 
eterize the recoil velocity: 



100 




VcM — 463 km s 



, hi) 



(9) 



If the recoil velocity of the coalescing binary exceeds 
the escape speed Wcsc, the holes will leave the galaxy 



^ We performed one calculation adopting the lFavata et all <2004') 
lower limits and checked that the fraction of displaced binaries 
amounts to about 5%. 



12 16 
redshift 



Fig. 3. — Averaged distribution of MBH mergers per merger 
tree of a lO-*^"^ Mq halo at ^ = 6. Solid line: when an accretion disc 
is present, the satellite MBH is dragged towards the central MBH 
on the viscous timescale [efficient merging). Dashed line: binaries 
evolve via DM-|-stars scattering {inefficient merging). 



Figure compares the MBH merger rates in two differ- 
ent models accounting for the orbital evolution of MBH 
binaries in the phase preceding emission of gravitational 
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waves. If gas processes (efficient merging) drive the 
MBH orbital decay, MBHs start merging efficiently at 
very early times, when host DM halos are still small. Al- 
though the absolute number of mergers is larger, more 
MBHs are ejected from DM halos due to the rocket ef- 
fect, because halos are smaller at earlier times. The two 
effects, increased mergers and increased ejections, com- 
pensate (see also figure O thus leading to similar SMBH 
masses in the main halo at z = 6. Note that the MBH 
merger rate in such a highly biased volume is not the 
average cosmic rate, as we are following only the hierar- 
chy leading to e xtrem ely high mass halos at z = 6 (cfr. 
iDiemand et alJ l|200!Tfl for similar considerations at lower 
redshift). MBHs in a high-overdensity volume experi- 
ence a much higher number of mergers with respect to 
the typical MBH evolving in an less overdense region. 
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Fig. 4. — Distribution of kicked BHs as a function the MBH bi- 
nary recoil velocity to the escape velocity from the host halo ratio. 
Solid line: fraction of MBHs with respect to the total number of 
seeds. Dotted line: only for binaries merging ai z < 10. Dashed 
line: number of MBHs over the number of binaries. The three pan- 
els refer to different models. Upper left panel: Eddington accretion 
rate only (e < 0.16), efficient merging. Upper right panel: Edding- 
ton accretion rate only (e < 0.16), inefficient merging. Lower left 
panel: MBHs are allowed to grow at supercritical rates, thus shift- 
ing the mass ratios of merging binaries towards mass ratios <C 1. 
Efficient merging is assumed. 



Figure^summarizes our results in terms of MBH bina- 
ries ejection. We have explored a large parameter space 
and we selected here a few representative examples. Ac- 
cretion and dynamical histories are strongly intermin- 
gled, as the kick velocity depends on the binary mass- 
ratio. The rocket recoil is substantial only for mass-ratios 
of order q > 0.1, which is believed to be the most proba- 
ble mass ratio for MBH binaries (due to the inefficiency of 
dynamical friction in bringing the satellite MBH within 
the infi uence of the main one in case of minor mergers, 
see also iSesana et al. (2005)). Interestingly, the Volon- 
teri & Rees model for supercritical accretion skews the 
mass ratios of MBH binaries towards lower values, as the 
conditions for supercritical accretion to happen are typ- 
ically fulfilled by only one of two merging systems. This 



latter result supports a biased and selective growth of 
high-z BHs, as the conditions for super-critical accretion 
appear to be fulfilled only in halos with Tvir > lO^K, 
which represent 3-(T or 4-cr peaks in the field of density 
fluctuations . If only a tiny fraction of the MBHs, those 
hosted in the most massive halos at this early time, un- 
dergo a rapid growth, the conditions for a mild effect of 
the recoil are naturally met, regardless of which detailed 
kick velocity calculation is considered. The distribution 
of merging MBHs mass ratios is in fact skewed towards 
low mass ratios { q — m2/mi <^ 1), where the expected 
kick velocity is low (Figure 



5.2. MBH accretion 

The mass growth history of a MBH ending up as a 
SMBH in a Mh = 10" Mq at z = 6 is shown in figureEl 
for a series of different models. MBH mass increase by 
mergers is taken into account, as well as mergers 'nega- 
tive feedback' onto the MBH growth due to the possibil- 
ity of binary ejection following coalescence. 




redshift 

Fig. 5. — Averaged mass growth history of the MBH at the 
center of a 10^^ Mq halo at z = 6. Upper left panel: Eddington 
accretion rate only (maximum radiative efficiency e = 0.16). Upper 
right panel: Eddington accretion rate only (maximum radiative 
efficiency e = 0.12). Lower left panel: super-critical accretion, in 
a fat disc with = 0.1. Lower right panel: Eddington accretion 
rate only (maximum radiative efficiency e = 0.16), and 25 times 
larger seed holes abundance. A cross marks the locus of SDSS 
1148-1-3251. 



As a reference, we have considered an upper limit to 
the radiative efficiency of e = 0.16, or d = 0.9 adopt- 
ing the standard conversion, inspired by iGammie et alJ 
(2004) simulations. Figure |21 shows, however, that MBH 
masses of order 10^ Mq can be reached by z = 6 only 
if dynamical effects are not too destructive. In fact, if 
we adopt the iBaker et all l)200(iD scaling (Eq. EJ for the 
gravitational recoil velocity, MBH masses are always less 
than a billion solar masses. The constraints imposed by 
z = 6 quasars can be met if either (i) the radiative ef- 
ficiency is lower than emax = 0.12 (corresponding to a 
spin parameter a — 0.8) or (ii) if the number of seed 
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MBHs is much more than 25 times^ larger (i.e. seeds in- 
habit density peaks ^ 3.5 — a at z=24. Note that in this 
case massive black holes would be expected also in dwarf 
galaxies with total masses well below Mh = 10" M© to- 
day). Alternatively, (iii) if the MBHs go through a phase 
of supercritical accretion. 

If the dynamical influence of the 'gravitational 
rocket' is bette r described by (e.g. iFavata et alJ 1200^ 
iDamour fc Cop akumar 2006), though, the observational 
constrains can be met more easily, allowing for higher ra- 
diative efficiencies or rarer seeds. It is worth noting that 
a radiative efficiency e > 0.16 can be accommodated only 
if the influence of mergers is large (see also Shapiro 2005 
and Yoo & Miralda-Escude 2004), e.g. seeds in lower 
(T-peaks and efficient merging of binary MBHs. 



(—2,-3.5). Letting the luminosity break vary between 



L*g^7 X 10"L 



o and L*g = 2 X lO^^L© leads to indis- 
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Fig. 6. — Luminosity function of 2 = 6 quasars, resulting from 
different models. The thick solid line shows the best fit LF from 
Fan et al. 2004. The dashed lines two model LF with faint and 
bright end slopes of the luminosity function (a, /3) = (—1,-2.5) 
and («, /3) = (-2,-3.5). Upper left panel: Eddington accretion 
rate only (maximum radiative efficiency e = 0.12). Upper right 
panel: Eddington accretion rate only (maximum radiative effi- 
ciency e = 0.16), Lower left panel: super-critical accretion, in a 
fat disc with = 0.1. 



We have compared theoretical luminosity functions 
(LF) with the most recen t determination of the quasar 
blue LF from the SDSS f Fan et all 1200^ . The SPSS 
sample s only the very bright end of the LF; iFan et aTl 
l)2001a() fit a single power-law to the data, which span 
luminosities larger than 10^'^ Lq in the blue band. Up- 
per limits on the faint end of the luminosity function, 
though, can be derived from the non-detec tion of 2: = 6 
quasa rs in the Chandra Deep Field-North ijBarger et al.l 
120031) and in the Ca nada-France High-z Quasar Survey 
(jWillott et alJl2005j) . We have a dopted here the pa,- 
rameterization of the LF given bv lWiUott et"al] 1I2OM) 
and compared our results with two representative cases 
for the faint (a) and bright {(3) end slopes of the lu- 
minosity function: (a,/3) = (—1,-2.5) and (a,/3) = 

Computational costs did not allow to explore initial conditions 
with a larger number of MBH seeds. 



tinguishable curves. Figure shows that quasars bright 
enough to reproduce the observed LF cannot be created 
under the assumption of accretion is limited at the Ed- 
dington rate, occurring via a disc which can spin up 
black holes up to a = 0.9 (corresponding using the stan- 
dard conversion to e = 0.16). The constraints from the 
faint end of LF are still very weak (cfr. Figure 5 in 
iWillott et a l. (2005j)), so we cannot rule out either a 
model with a lower Cmax or a model with supercritical 
accretion based on these results only. 

6. DISCUSSION 

The detection of luminous quasars at z ~ 6, suggest- 
ing an early growth of SMBHs, deserves a special inves- 
tigation within the hierarchical scenario of galaxy forma- 
tion. Following a series of papers tracing the seeds of the 
SMBH population observed today from the first stars, we 
have here focused on the constraints set by the presence 
of SMBHs already in place when the Universe was less 
than one billion year old. 

The haloes that we choose for our investigation, Mh > 
10^^ Mq at z = 6, are 5.5-cr density fluctuation. They 
are not therefore representative of the typical halo mass 
at z = 6. On the other hand, the masses are chosen 
requiring that their number density matches the number 
density of quasars derived by Fan et al. (2004). That 
the highest massive halos experienced strong evolution at 
early times is in line with the 'anti-hierarchical' evolution 
of galactic s tructures sugges ted by high-redshift galactic 
surveys (e.g. lKodama et al] l2004). A simple exercise can 
help us quantify the global importance of the black holes 
hosted in the selected halos. 

If the correlation between black hole mass and halo 
circular velocity holds up to z = 6, we can use the Press 
& Schechter formalism to determine the MBH density 
as a function of halo mass at z = 6. The MBH mass 
density in halos with mass larger then Mh — 10^^ Mq 
is ~ 20MqMpc"^, while for Mh = 10^^ Mq we have 
~ 10^MqMpc~3, and for Mh = 10" M© we have ~ 
lO^MoMpc-^. So, if all halos with, e.g., Mh = 10" Mq, 
at z = 6 host a SMBH whose mass scales with the local 
-^•^BH — Vc , the SMBH density at z = 6 is similar to that 
at z = 3, leaving very small room for accretion between 
z = 6 and z = 3. This implies that either the A/bh — Vc 
relation is redshift dep endent, and in the o pposite sense 
than suggested by the I Walter et al.l lj2004[l observation, 
or the occupation fraction of black holes is not unity at all 
redshift. We define the MBH occupation fraction, as the 
average number of halos hosting a MBH over the total 
number of halos in a co moving volume. Volonteri, Haardt 
& Madau 2003 (see also lMarulli et all l|2006) ) show that 
models in which SMBHs evolve from seeds forming in 
biased regions (e.g., remnants of PopIII stars), predict 
a SMBH occupation fraction which decreases with halo 
mass. 

If the occupation fraction is a declining function of the 
halo mass, we expect MBH-MBH interactions to be more 
common for massive halos. MBHs hosted in low mass 
halos grow in isolation, while MBHs lurking in massive 
halos have a much higher probability of experiencing a 
MBH-MBH merger. So, on one hand MBHs hosted in 
low mass halos have a larger probability of being ejected 
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by a kick, on the other hand their probabihty of ex- 
periencing a merger in the first place is smaUer. We 
might therefore expect that a SMBH can grow almost 
undisturbed by ^ = 6 up to a billion solar masses, if 
hosted in a, say, 10'^"'^ M0 halo. A requirement in this 
case is, however, that the mass of the black holes, in 
these high redshift systems, can grow to masses much 
larger than allowed by the local Mbh — Vc or correlation. 
This i s the opposite of wha t seen in numerical simula- 
tions (Robertson ct al. 2005), but in agreement with the 
[Walter ct al. ( 2004) results. 

Summarizing, if the SMBHs powering z = 6 quasars 
are hosted in halos with Mh ^ 10^^ M©, they can grow 
almost in isolation, and accretion can be limited to the 
Eddington rate, provided the radiative efficiency is not 
too high, but the mass of the SMBHs must be allowed 
to exceed the nominal Ms/f— velocity dispersion correla- 
tion. 

If the host halo masses are of order Mji 2± 10^'^ Mq, 
dynamical interactions are important in their evolution- 
ary history. In this case, either kick velocities are at the 
lower end of the predicted range, or accretion occurs at 
supercritical rates for some periods of time. 

If, instead, high-redshift MBHs can accrete at super- 
critical rate during an early phase, reproducing the ob- 
served MBH mass values is not an issue, even in the 
case that the recoil velocity is in the 'upper limits' range. 
Similar conclusion can be reached also with alternative 
models for early supercritical accretion onto MBH seeds 
formed out of direct collapse of gas in high-redshift halos 



(Begelman, Volonteri & Rees, 2006). 

The emerging picture, therefore, points to a Darwinian 
natural selection scenario. At early times only a small 
fraction of MBHs, those hosted in the highest density 
pea ks (Tyir > lO^K), ev olve rapidly and efficiently (see 
also lJohnson fc Bromml QOOG)). The effects are twofold: 
not only these rare BHs become quickly supermassive 
owing to efficient accretion, but also the effects of danger- 
ous dynamical interactions are softened, as the merging 
of low mass ratio MBH binaries is favored. Alte rnative 
mode ls for MBH seed formation (e.g. .Koushiappas et alJ 
120041) tend to select as sites of black hole formation halos 
with Tvir > lO^K as well. 

Clarification of the early evolution of SMBHs can come 
from detection of the gravitational wave signal from their 
inspiral, as different models for the MBH early evolution 
predict different mass-ratio distributions for the merging 
MBHs. Reliable calculations narrowing the uncertain- 
ties on the kick velocity due to gravitational emission for 
spinning MBH binaries will help clarify the threat that 
kicks represent for the early evolution of the MBH pop- 
ulation. 



We have benefitted from useful discussions with Picro 
Madau, Francesco Haardt, Cole Miller, and Monica 
Colpi. MV thanks the Department of Astronomy & As- 
trophysics of the University of California at Santa Cruz 
for the warm hospitality. We would like to thank the 
referee for helpful comments. 
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